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Abstract 



We study a class of tridiagonal matrix models, the "q-roots of unity" models, 
which includes the sign (q = 2) and the clock (q = oo) models by Feinberg and 
Zee. We find that the eigenvalue densities are bounded by and have the symme- 
tries of the regular polygon with 2q sides, in the complex plane. Furthermore 
the averaged traces of M'^ are integers that count closed random walks on the 
line, such that each site is visited a number of times multiple of q. We obtain 
an explicit evaluation for them. 
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1 Introduction 



Random matrix ensembles are extensively studied since the early works of 
Wigncr and Dyson, as effective models for the description of statistical proper- 
ties of spectra of complex physical systems, which include resonances of heavy 
nuclei, quantum billiards, mesoscopic transport and quenched QCD - Q. 

Tridiagonal matrices with random entries naturally occur in the simplest 
models of disordered one-dimensional crystals, beginning with the works by 
Dyson and Schmidt If the spring constants are fixed while the masses 
are random, one is led to a tridiagonal Hermitian matrix with random entries 
on the main diagonal, often referred to as a random site problem. 
More recently, several authors studied the density of eigenvalues of non-Hermitian 
tridiagonal matrices still having the random entries in the main diagonal, af- 
ter the model introduced by Hatano and Nelson to describe the motion of 
vortices pinned by columnar defects in a superconductor . In this model, the 
spectral density was obtained only in the case of a Cauchy probability distri- 
bution for the independent random entries js) , : in the large- iV limit, the 
eigenvalues lie on a " squeezed ellipse" in the complex plane with the addition of 
two "wings" on the real axis, which appear provided the strength of the random 
site entries exceeds a critical value . It was known for a long time ||T^ that 
complex non-Hermitian random matrices may have eigenvalues filling a two- 
dimensional area in the complex plane. For a class of tridiagonal non-Hermitian 
models, with random hopping and random sites, Goldsheid and Khoruzhenko 

found conditions for the probability distributions, such that in the large- A'^ 
limit the eigenvalues converge to a curve in the complex plane. Yet, in more 
general settings it seems still impossible to predict whether a non-Hermitian 
ensemble of random matrices has eigenvalues converging to a curve or filling 
a two-dimensional area or a fractal. These difficulties suggest the usefulness 
of investigating different types of tridiagonal non-Hermitian random matrices, 
possibly by a variety of techniques, even for models not directly related to a 
physics problem. Feinberg and Zee fl^ considered a class of non-Hermitian 
random hopping models, with eigenvalue equation 

tj-iipj-i + s*'ipj+i = Eipj (1.1) 

and given probability distributions for the hopping amplitudes tj , Sj. In their 
first model, called the "clock model" , the hopping amplitudes are independent 
random phases. They found , numerically, that the eigenvalues are distributed 
in a disk in the complex plane, centered at the origin, with rotational invariance. 
In a second model , called the "sign model" , hopping amplitudes are indepen- 
dent and randomly equal to ±1. In the case of open chain, it is easy to see that 
the eigenvalues are functions of the — 1 products tj s* , therefore , for the open 
chain, the distribution of eigenvalues of the " sign model" is also obtained by the 



study of eigenvalues of the matrix M{x,l) , given below in eq.(2.1), where the 
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N — I independent variables independent and randomly equal to ±1. 

In this paper we study a more general model , the "q-roots of unity" model, 
where each independent variable Xi is one of the q roots of unity, with uniform 
probability. Of course, for g = 2 , it is again the "sign model". 
For any q , the sum of the moduli of the two non-zero entries in any row of the 
matrix M{x, 1) equals two, therefore by Gershgorin circle theorem, all eigen- 
values are in the disk \E\ < 2. We find that for any N the eigenvalues of the 
"q-roots of unity" model are inside the regular polygon with 2q sides. The 
derivation of this very unusual boundary, together with the symmetry proper- 
ties of the eigenvalue distribution are of intrinsic interest and are given in Sect. 2. 
In the limit q oo , which corresponds to assuming the random variables xt 
to be random phases, the support is the disk centered in the origin, with radius 
equal two and the density only depends on \E\. This case is the "clock model" 
of ref. Our second result is the evaluation of the moments < tr Af" > by 
mapping it into the problem of counting the returns to the origin of a restricted 
class of one-dimensional random walks, in Sect. 3 jl^ . The number of random 
walks in one dimension that originate at a given point, which we call the origin, 
and after 2n random steps of unit length to the right or to the left, return to 
the origin (not necessarily for the first time) is (2n)!/(n!)^. However, if we select 
among them the walks where each site different from the origin is visited an even 
number of times, the walks have to consist of a number of steps multiple of 4 
and their number is smaller. Let us call such random walks the even- visiting 
walks. In the limit ^ oo , the number of even-visiting walks provides the 
moments < tr Af " > of the " sign model" . In the same way, the number of 
returns to the origin of one dimensional walks where each site of the walk is 
visited a number of times multiple of q provides the moments < tr Af" > of 
the "q-roots of unity" model. 

As we mentioned , we do not know physics models related to the "sign 
model" . It may then be proper to spend a few words to justify the present 
investigation. The evaluation of moments < tr M" >, where the tridiagonal 
matrix M{x, 1) is given in eq. (|2.l[) , is conveniently mapped into a random walk 
problem in one dimension, independently of the probability densities of the ran- 
dom variables x^. Furthermore, provided the Xi are real independent variables 
and the probability densities are even functions, only the even- visiting walks are 
relevant. However, only for the "sign model" the evaluation of the moments < 
tr Af" > corresponds to counting the even- visiting walks, whereas if a different 
even probability density is considered, to each even-visiting walk one has to 
assign a weigth. 
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Our main result , for the "sign model" is 

k 

lun 1 < trM^^- E E Sln,.n,,-,n.] (1.2) 

where the sum X]{n } over the t positive integers rii with the restriction 
rii + n2 + ■ ■ ■ + fit = k and S'[„^.„,.....„j] is the number of the relevant walks with 
"width" w^t : 

A very similar result, for arbitrary value of q, for the "q-roots of unity" model, 
is obtained in Sect. 3. 

As this work was completed, a new work appeared [jisl where the spectral 
density of the "clock model" was obtained. Since our work does not provide 
quantitative information on spectral densities, it has little overlap with it. How- 
ever several very interesting assertions , like the analysis of the two non-linear 
maps r+ , r_, and the existence of bound states , discussed for the "clock 
model" , are equally valid for the " sign model" . 



2 The q-roots of unity model 

We consider the ensemble £n {q, 1) of tridiagonal random matrices of size NxN, 



M{x,l) = 
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(2.1) 



where the random entries Xi, i = l,2,..,iV — 1, are g— roots of unity: xj = 
1. More precisely, they are independent and identically distributed random 
variables, with uniform probability distribution over the set of the g— roots of 
unity: 



P{x) 



n ^ — ^ 



(2.2) 



The ensemble consists of q^ ^ different matrices. 

For later convenience, we introduce the notation M{x,y) to denote a general 
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tridiagonal matrix M of order with the upper and lower diagonals specified 
by vectors x = {xi, ..,xn-i} and y = {yi, ..,yN-i}. 

In this section we are interested in the investigation of the symmetries and 
boundaries of the eigenvalue distribution of the ensemble, in the large N limit. 
We consider , for finite N, the set cr(g, TV) in the complex plane of all the eigen- 
values of matrices belonging to the ensemble and show that it has nice symmetry 
properties. 

Proposition 1 : The set a{q, N) is invariant under the transformations 
i) -E, ii) E^E*, iii) E Ee'^ 
Proof: the characteristic polynomial of a generic tridiagonal matrix M{x, y) only 
depends on the products xiyi, i = 1,2, .., N — 1. Therefore, for a given matrix 
M{x, 1) in the ensemble we introduce the two matrices M{z, z) and M{—z, —z) 
where zf = Xi. Clearly they have opposite eigenvalues, meaning that both E 
and —E are eigenvalues of M(x, 1). 

The ensemble £N{q, 1) is closed under complex conjugation, and if {Ei} are the 

eigenvalues of M{x, 1), {E*} are the eigenvalues of M{x* , 1). 
The set of matrices M(z, z), where each Zi is randomly chosen in the set of the 
2g— roots of unity, is invariant under multiplication of a matrix by the scalar 
giTT/g Therefore, if {Ei} are the eigenvalues of M(a;, 1), {£^,6*'^/'} are the eigen- 
values of M{xWq, 1). 

We conclude that the set a{q, N) has the above stated symmetries (note that 
invariance i follows from iii). • 

Proposition 2 : The set (T{q, N) is contained in the regular 2g— polygon in 
the complex i?— plane, with corners Ek = 2e*'^'^/'?, k = 0, 1, .., 2g — 1. 
Proof. Let E be an eigenvalue of a given matrix M{x, 1) in the ensemble; 
because of the symmetry, we assume that < aigE < ir/q. The eigenvalue 
equation for M{x, 1) is Ui-i + XiUi+\ = Eui, with uo = un+i = 0. Let Ik be 
the set of integers j = 1, .., — 1 such that xj = w'^. After multiplication by 
u* of the eigenvalue equation and summation over i: 

N N 

E = u*Ui-i + ^ XiU*Ui+i = 

i=l i=l 

JV q 

i=l k=l iQlk 

q 

= XI XI "i+i"* + wlu*Ui+i = 
fc=i ie/fc 

q 

fe=i ielk 
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We finally take the real and imaginary parts of E and build the majorization: 

q 

ReEcosf— ) +Im£'sin(— ) = cos(k- - — ) (<Wj_ie~'''? +c.c. ) 

W '2q' l^^ ^ q 2q' {^\^ J 

fe=i ^ ^ ie/fe ^ fe=iie/fe 

< 2C0S(;^) 

Therefore: ImE < cotg(^)(2 — BjeE), which ends the proof. • 

Prom the discussion of proposition 1, we conclude that the distribution of 

the eigenvalues p{E,E*) for the ensemble £i^{q,l) has the listed symmetries, 
and therefore it is a symmetric function of the variables -E^"^ and (E*)'^'^. This 
puts contraints on the moments. 

Proposition 3 : {E"^E*^) ^ only if m + n is even and \m — n\ = 2rq. 
Proof: Since the expectation value is symmetric in m and n, we put m > n. 
Due to the symmetry E — > ^E of the density, non- vanishing of the expectation 
value requires m + n = 2s, so that m — n = 2£. By the 2g— fold symmetry of 
the density and the support a, we may evaluate the expectation by integrating 
over u), the angular sector < argE' <iT/qoi a: 



2q-l 



* ^—iklT / q\—£ 



{\Ef'E'^{E*)-'^) = j <fEp{E,E*)\E\^'^{Ee'''''/'iY{E 
•'•^ fe=o 

Wo get a non-zero result for the sum only if ^ = rq. • 

In the limit g' ^ oo, where each number Xi is an arbitrary number on the 
unit circle, the set cr(oo, N) is invariant under complex rotations, and the eigen- 
value density is a function of 



3 The moments and counting the returns of the 
relevant walks. 

In this section we study a set of moments of the probability density p{E,E*) 
of the "q-roots of unity" model, which have the interesting interpretation as 
counting numbers of the random walks on the line that return to the origin and 
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Figure 1: One of the even- visiting random walks, returning to site r after 8 
steps, with width w = 2. 



visit each intermediate site a number of times which is multiple of 2q. These 
moments have generating function 

where the expectation value is evaluated on the independent identically dis- 
tributed random variables Xi : 

N 

< f{M) >= Yl P{xi) dxi f {M{x, 1)) (3.2) 
The formal perturbative expansion of G{z) is 

k=0 

Since we are interested in the limit N —>■ oo, any term on the diagonal 
< [M'']r-r > has the same value, and the trace merely cancels the 1/A'' factor. 

It is useful to consider the one to one map between the non vanishing terms 
contributing to [M'^]a(, and the random walks in one dimension which originate 
at site a to arrive at site b after k steps. 

Let us consider the term fc = 8 

^ < MraM^bMbcM^dMdeMefMfgMgr > (3.4) 

a,byC^d,e,f,g 



7 



01 2345678 



Figure 2: One of the random walks not belonging to the class of even-visiting 
random walks. 



By recalling that the non zero matrix elements are Mij = 1 if j i — 1, and 
Mij = if j = i + 1 , each term in the sum (3.4) corresponds to a walk of 8 



steps, originating and ending at site r, with 4 steps up and 4 steps down. For 
instance, the sequence 

Xr+1 • 1 • 1 - Xr • 1 = 1 is shown in Fig.l, while the sequence Xr-Xr+i 

'1*1* Xf ' 1 * Xf ' 1 

shown in Fig. 2 . Each walk corresponding to a product of random 
variables rij(^j)"^ where all the powers nj are multiple of q yields a contribution 
+ 1, while the walks where at least one power nj is not a multiple of q are 
averaged to zero , therefore being a class of irrelevant walks. Then jj < tiM'' > 
is equal to the number of relevant walks of k steps, that is the walks from a fixed 
site r to the same site r , such that each intermediate site is visited a multiple 
of q times. 

Because the number of steps up (each associated to a random variable Xi) 
is equal to number of steps down (each associated to a factor one) , the total 
number of steps k of the relevant walks is a multiple of 2q and we rewrite eq.( |3.3[) 
as 

oo 

G(z) = y_^ ; Cfc = lim — <trM2«'= > (3.5) 

k=0 

This power series is absolutely convergent for \z\ > 2. In the case of g = cx) we 
obtain 

G{z) = ^ , for \z\ > 2 (3.6) 



Next we present the combinatorial evaluation of the coefhcients c^, for the case 




Figure 3: Here are shown the 14 paths of 8 steps, relevant for the "sign model" 



q = 2 , that is the number of returns to the origin for even visiting walks. The 
case of general integer values of q requires only trivial generalization, presented 
in the following paragraph. 

Counting the returns to the origin of even visiting walks . We shall indicate 
with N{2n-s, 2n_s+i, • • • , 2n_i ; 2no, 2ni, . . . , 2nj) the number of even visiting 
walks from site r to site r corresponding to the product 

(a;^-^)^"- • • • {xr-i)^"-^{xrf"'>{xr+i)^"^ ■ ■ ■ (xr+j)^"' ■ Ah the integers , ...rij 
are positive ; the semicolon separates the exponents rn associated to sites i < r 
from those associated to sites i > r. We omit strings of zeros external to the 
string of positive integers, just keeping one zero for walks visiting only sites 
s > r or only sites s < r . In the first case we write the multiplicity as 
A^(0 ; 2nQ, ■ ■ ■ ,2nj) , in the latter case wc write A^(2n_s, ■ • • , 2n_i ; 0). It is 
easiest to begin with the evaluation of iV(0 ; 2no, • ■ • , 2nj) . The length of the 



9 





































/ 


\ 




























/ 





























Figure 4: The figure illustrates the insertion of a 4-steps part over anyone of the 
tops of a previous relevant walk of width w, thus obtaining new relevant walks, 
of width li; + 1 , according to eq.(B.7) 



walk is 4(no + rti + • • • + The "maximum site" visited is the site r + j + 1 , 
visited luj times ; the "minimum site" visited is the site r , visited 2rto + l times. 

The number 7V(0 ; 27io, 2rti, • • • , 2nj, 2njj^\) is related to iV(0 ; 2no, 2ni, • • • , 2nj) 
in the following way : new walks of length two corresponding to (xr+j+i • 1) may 
be inserted in each of the maxima of the previous walk. Since 2nj+i identical 

objects are placed in 2nj places in ways, we obtain 

iV(0 ; 2no, 2ni, • • • , 2nj,2nj+i) = 
= f^''-'%+^''-'~MA^(0;2no,2ni,...,2n,) (3.7) 
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By iterating eq.(3.7) with the initial condition iV(0 ; 2no) = 1 one obtains 



A^(0; 2no,2ni,- 



,2n. 



(=0 



2ni+i + 2ni 



1 



(3.8) 



We proceed to evaluate the multiplicity of the relevant random walks that 
visit sites r, r + 1, • • • , r + j as often as before, and in addition visit 2n_i times 
the site r — 1 . Each walk of this class may be obtained by inserting 27i_i walks 
of length two (1 • Xr-i) in each of the 2no + 1 minima of the walk of the previous 
class. Therefore 



iV(2n_i; 2nn,2ni, - ■ ■ ,2nj) 
2n_i + 2nQ 
2ri_i 



N{0;2no,2ni,...,2nj) (3.9) 



The procedure may be repeated to include the relevant walks which visit sites 
r — 2, r — 3, • • • , r — s. We obtain 



iV(2n_s, 2n^s+i, ■■■ , 2n_i ; 2no, 2ni, • • • , 2nj) 

s-2 



-j-j- j 2n^s+p 

.p=o 



2n_ 



S+P+ 



1-1 



2n_ 



s+p 



2n_i + 2no 
2n_i 



n 



2n,;+i + 2n,; — 1 



2n 



i+l 



(3.10) 



The coefficient Ck is the number of the even visiting walks of 4fc steps and 
it is the sum of the multiplicities N{2n^s, 2n-^s+i, • • • , 2n_i ; 2no, 2ni, . . . , 2nj) 
given above, where k — n^s + '^--s+i + ■ ' ' + "-j- 

The evaluation may be someway simplified, by considering walks of fixed 
width, that is the difference between the "maximum site" visited and the "min- 
imum site" visited . We consider the set of ordered partitions of k into positive 
integers [ni, ^2, • • • , nt] where fc = "-p- Each ordered sequence , [rii, ri2, • • • , nt] 
, corresponds to t + 1 classes of even visiting walks, which are associated to the 
products 

(:E,_,)'"H^r-t+i)'"^---(^r-i)'"* ; 

(x,_t+i)2"H^r.-t+2)'"^ • • • (a^r)'"' ; 



(x,)2«i(a;,+i)2«^...(a;,+,_i)2"* 



(3.11) 



All walks in eq.(3.11) have the same width w = t. Their multiplicities, given in 
eq.(3.10) are simply related and their sum is 



S[ni,n2,---,nt] ~ „ H 
"1 ._i 



2k W ^ 2n,+i + 2n.i - 1 



2n. 



■i+l 



S[ni] — 2 



(3.12) 
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Next we sum over the ordered partitions [tii, 712, • • • , n*] of fc into t parts and 
finally over the difi^erent widths, from 1 to fc 



Ck 



H ^l" 



(3.13) 



where the sum X]{ri } over the t positive integers with the restriction 
ni + n2 + ■ ■ ■ + nt = k. The evaluation of eq.( 3.13| ) may be automated and we 
find the first coefficients Ck '■ 



Co 


1 


Cl 


2 


C2 


14 


C3 


116 


C4 


1 110 


C5 


11 372 


C6 


123 020 


C7 


1 384 168 


C8 


16 058 982 


Cg 


190 948 796 


ClO 


2 317 085 924 


Cll 


28 602 719 576 


C12 


358 298 116 092 


Cl3 


4 545 807 497 272 


Cl4 


58 321 701 832 408 


Cl5 


755 700 271 652 816 


C16 


9 878 971 460 641 414 



(3.14) 



The ratios c„/cn_i rise monotonically with a rate slower at higher values 
of n. We know from the previous section that the eigenvalues of the random 
matrix M are inside the square with vertices in the four points ±2 , ±2i , it 
then follows that c„ ~ a" , with a < 16 . The plot of the coefficients c„ for 
n = 4 to 19 in Fig. 5 is consistent with the expected asymptotics. • 

Counting the returns to the origin of walks for generic q . The combinato- 
rial evaluation given above holds with trivial modifications for the "q- roots of 
unity" model, described in eqs.(2.1) , (2.2). Again tr M'^ corresponds to one- 
dimensional walks returning to the original site after k steps, each walk being 
the product of fc/2 random variable s Xj . The average of each product, with the 
factorized probability distribution (|2.2| ) vanishes unless each random variable 
occurs with a power multiple of q , in this case the product has average one. 
Therefore fc/2 is a multiple of q and the class of relevant walks is such that each 
site (apart the origin of the walk) is visited a number of times multiple of q. 
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Figure 5: The triangles are the values of (l/4fc) log2 Ck versus (l/4fc), for k — 4 
to 19. The convergence for large values of k is consistent with ^ ^ with 
a < 16 



Let A'^(0 ; qnQ, qni, qrij) be the number of the relevant walks corresponding 
to the multiphcity of the product (x^)*"" (xr+i)'"^ . . . (xr+j)*"'. The insertion 
above the top sites or below the bottom ones proceeds as before leading to a 
trivial generalization for the number of relevant walks S[ni,n2,---,nt] having fixed 
width t 

2fc / qui+i +qn,-l\ 



k 

Cfc = X! X! "^["i. "2, ■■■,«*] (3.16) 

where the sum X){n } over the t positive integers Ui with the restriction 
ni + n2 + ■ ■ ■ + nt — k . • 

Let us close with few words for the special value q = 1. In this case the 
matrix Al in eq.( |2.lD is no longer random, it is the real symmetric matrix with 
all Xi = 1. Its spectral distribution is well known and the analysis in random 
walks is not necessary. Still one may find that — tr Af'^(l,l) = i^k] is 
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the number of one-dimensional random walks returning to the origin after 2k 
steps. We checked , for several values of k, that eqs.( |3.15 ) , (3.16) reproduce 
the correct value for q = 1. 
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